Magnetovolume effect in Ce(Nii_^Cu2;)5 alloys 
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Magnetic susceptibility x the isostructural Ce(Nii_a;Cu2,)5 alloys (0 < a; < 0.9) was studied 
as a function of the hydrostatic pressure up to 2 kbar at fixed temperatures 77.3 and 300 K, using 
a pendulum-type magnetometer. A pronounced magnitude of the pressure effect is found to be 
negative in sign and strongly (non-monotonously) dependent on the Cu content, showing a sharp 
maximum at a; ~ 0.4. The experimental results are discussed in terms of the valence instability 
of Ce ion in the studied alloys. For the reference CeNis compound the main contributions to x 
and their volume dependence are calculated ab initio within the local spin density approximation 
(LSD A), and appeared to be in close agreement with experimental data. 

PACS numbers: 71.20.Eh, 75.30.Mb, 75.80.+q 
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I. INTRODUCTION 

Many of Ce intermetallics are characterized by a strong 
hybridization of the magnetic 4/ electrons with the con- 
duction electron states resulted in delocalization of the 
4/ level and a change of its occupancy, and hence the 
Ce valence. As is evident from measurements of X- 
ray absorption and lattice parameters [l|, together with 
the magnetic 0, [1], electric and thermoelectric proper- 
ties [3|] in the isostructural Ce(Nii_a.Cua;)5 alloys, the Ce 
valence decreases steadily from Ce^"'' to Ce'^+ with in- 
crease of the Cu content, and the system undergoes a 
series of transitions from the nonmagnetic metal with 
empty 4/ level {x — 0) through the intermediate va- 
lence (IV) state combined with a nonmagnetic dense 
Kondo state (0.1 < a; < 0.8) to the magnetic 4/ metal 
(0.9 < 2^ < 1). Thus, the reference CeNis coinpound 
is the exchange-enhanced itinerant paramagnet [l|, 0) 0] 
with the temperature dependent magnetic susceptibility 
exhibiting a broad maximum around 100 K, similar to 
those observed in YNis and LuNis 0, 0] ■ On the other 
side, the CeCus compound is a Kondo lattice antifer- 
romagnet with Tn = 3.9 K and Tk = 2.2 K ?|. The 
magnetic susceptibility in CeCus at T > 50 K obeys 
a Curie- Weiss law with the effective magnetic moment 
value close to that expected for Ce'^^ state 0, 0] ■ Due 
to a direct relation between magnetic properties and the 
rare earth (RE) valence state, and also the strong correla- 
tion between the valence itself and RE ionic volume, the 
RE compounds with unstable / shell exhibit a large mag- 
netovolume effect. Therefore, a study of pressure effect 
on magnetic properties of the systems with variable RE 
valence is of great interest to gain insight into a nature 
of the IV state. 

Here we report results of our investigation of 
the pressure effect on the magnetic susceptibility of 



Ce(Nii_a;Cua;)5 alloys in a wide range of Cu concentra- 
tions. The experimental results are supplemented by cal- 
culations of the magnetovolume effect value for the refer- 
ence CeNis compound, using a modified relativistic full 
potential approach within linearized " muffin-tin" orbital 
method (FP-LMTO). 



II. EXPERIMENTAL DETAILS AND RESULTS 

The polycrystalline samples of Ce(Nii_2;Cu2;)s alloys 
(0 < X < 0.9) were prepared by arc-melting of a stoi- 
chiometric amount of initial elements in a water cooled 
crucible under protective argon atmosphere. The study 
of X-ray powder diffraction at room temperature re- 
vealed that all samples crystallize in CaCus-type hexag- 
onal structure, and obtained data on their lattice param- 
eters agree closely with that published in literature. Any 
other phases were not detected within the resolution of 
the X-ray technique. 

The pressure effect on the magnetic susceptibility x 
was measured under helium gas pressure up to 2 kbar at 
two fixed temperatures, 77.3 and 300 K, using a pendu- 



TABLE I: The magnetic susceptibilities and their pressure 
derivative for Ce(Nii_i,Cui,)5 alloys at 77.3 and 300 K. 
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FIG. 1: Pressure dependence of the magnetic susceptibility 
of Ce(Nii_^Cu:„)5 alloys at T=77.3 K (a) and 300 K (b) nor- 
malized to its value at P = 0. 



lum magnetometer placed into the nonmagnetic pressure 
cell The relative errors of our measurements, per- 

formed in the magnetic field H = 1.7 T, did not exceed 
0.05%. 

In Fig. 1 the typical pressure dependencies of the 
magnetic susceptibility for Ce(Nii_2;Cua;)5 alloys demon- 
strate a magnitude of the pressure effect ant its linear 
behavior. For each temperature the values of x s-t am- 
bient pressure and their pressure derivatives, dlnx/dP, 
are listed in Table I. Both values include corrections for 
a weak field dependence of x caused by ferromagnetic 
impurities, which are less than 5%. The negative sign 
of the pressure effect is consistent with anticipation that 
high pressure has to increase the valence, since the Ce ion 
in the higher valence (less magnetic) state has a smaller 
volume. 

Of particular interest is a strong and non-monotonous 
concentration dependence of the pressure effect which 
shows a sharp maximum in vicinity of x ~ 0.4 for both 
temperatures, 77.3 and 300 K (Fig. 2(a)). A comparison 
between the obtained experimental results and the data 
on concentration dependence of the lattice parameter a 
and the effective Ce valence from Ref. [H (Fig. 2(b)) 
indicates that the maximum in dlnx{x , T) / dP corre- 
lates with a drastic change of a (and i>) around x ~ 0.4 
(i/ ci^ 3.5). 

It is interesting to note that a similar peculiarity in 
dlnx/dP versus valence was observed for various Yb 
compounds at room temperature [llj . As was shown, 
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FIG. 2; (a) Pressure derivative of the magnetic susceptibil- 
ity dlnx/dP in Ce(Nii_^Cux)5 alloys at 77.3 and 300 K. (b) 
Deviation of the a lattice parameter, Aa, in Ce(Nii_2,Cu2,)5 
alloys from the a(x) dependence for the Ce ion assumed to 
be in a trivalent state (left scale) and the Ce valence deduced 
from X-ray absorption studies (right scale) at room ternpera- 
ture versus Cu content x (according to the data of RefU). 



the relative change of x with pressure is the most pro- 
nounced also at the half-integer value of valence, v ~ 2.5, 
but contrary to the Ce compounds, it has a positive sign, 
as can be expected. 



III. THEORY FOR CENI5 



Ah initio calculations of the electronic structure were 
carried out for the reference compound CeNis by em- 
ploying a modified FP-LMTO method [H, The 
exchange-correlation potential was treated in the LSDA 
approximation [l^ of the density functional theory. To 
analyze the observed magnetovolume effect value in 
CeNis, the magnetic susceptibility and its volume de- 
pendence were calculated within the modified method, 
wherein the external magnetic field H was taken into ac- 
count by means of the Zeeman operator, i7(2s-fl). The 
latter was incorporated in FP-LMTO Hamiltonian 
for calculations of the field-induced spin and orbital mag- 
netic moments. The corresponding contributions to mag- 
netic susceptibility were derived from the field-induced 
moments, which have been calculated in an external mag- 
netic field of 10 T. 

The electronic structure calculations were performed 
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for a number of lattice parameters close to the experimen- 
tal one (the ratio c/a was fixed at its experimental value). 
The equilibrium lattice spacing ath and corresponding 
theoretical bulk modulus Bth were determined from de- 
pendence of the total energy on the unit cell volume, 
E{V), by using the well known Murnaghan equation, and 
appeared to be Oth = 8.96 a.u. and Bth = 1-9 Mbar. The 
Murnaghan equation is based on the assumption that 
the pressure derivative B' of the bulk modulus B is con- 
stant. By using the evaluated from the Murnaghan equa- 
tion value of B' = 3.73, we have estimated B^^ — 1.45, 
corresponding to the experimental Oexp = 9.2 a.u. 
This estimation appeared to be in nice agreement with 



the available experimental value. Be 



1.43 Mbar hf 



The differences between the equilibrium theoretical ath 
and Bth and experiment on bulk properties of CeNis are 
presumably related to the overbonding tendency of the 
LSDA approach [l2|. 

The strongly volume dependent spin contribution to x 
originates predominantly from the 3d states of Ni. Re- 
garding the orbital contribution to X: it comes mainly 
from electrons in the atomic sphere of Ce and amounts 
to about 20% of total susceptibility. At the theoreti- 
cal lattice parameter, the calculated total susceptibility 
(2.9 X 10~^ emu/mole) appeared to be very close to the 
experimental value (3.0 x 10~^ emu/mole at T = K 
The calculated volume derivative of susceptibil- 
ity, dlnx/dlny = 4.2, is in agreement with that resulted 
from the experimentally observed pressure derivative for 
CeNis at T = 77.3 K, dlnx/dlnl/ = 3.9 ± 0.4. Thus it 
has been demonstrated, that LSDA provides an adequate 
description of the strongly exchange enhanced magnetic 
susceptibility of CeNis and its pressure dependence. 



IV. DISCUSSION 

As is shown, the LSDA allows to describe the mag- 
netovolume effect in the reference CeNis compound that 
gives grounds for future application of ab initio LSDA ap- 
proaches to some Ce(Nii_a;Cua;)s alloys. Here, however, 
we shall restrict our consideration of the experimental 
data in alloys within a phenomcnological approach. 



A. Concentration dependence 

Anticipating the pressure effect on the magnetic sus- 
ceptibility to arise mainly from the change of Ce valence 
i/, or the fractional occupation of the 4/^ magnetic state 
ri4f {v — A — n4f), the pressure effect can be analyzed 
within a simple relation 



dlnx(r) dlnxiT) dn^t 



dP 



dn. 



■4f 



dP 



(1) 



in terms of the pressure dependence of n4f (or u). The 
most reliable results of such analysis would be expected in 
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FIG. 3: (a) Values of Inx at 77.3 K. (b) dna/dP plotted 
against n^f for Ce(Nii_a;Cua;)5 alloys, o and A denote the 
data obtained within Eq. ((T} and Eq. @, respectively. Points 
for X — 0.45 are interpolation of the experimental data on 
concentration dependence of x and dlnx/dP. 



the Cu-rich alloys at low temperatures where the 4/ con- 
tribution X4f becomes dominant (x ~ X4f)- In Fig- 3(a) 
the X versus n^i dependence is shown for Ce(Nii_a;Cua;)5 
alloys (0.4 < a; < 1) at 77.3 K, which was obtained by 
using the experimental xi^) values from Table 1 and 
iy{x) data of Fig. 2(b). A substitution of the resulted 
from Fig. 3(a) derivatives 51nx/9«4f and experimental 
data on dlnx/dP at 77.3 K into Eq. ^ gives the value 
dn4f/dP which strongly depends on i/4f (Fig. 3(b)). As is 
seen, the maximum value of dn4f is expected at v^f = 0.5 
(j^ — 3.5) to be about —6.5 ± 1.5 Mbar~^. The corre- 
sponding estimates of the valence change under pressure, 
di^/dP = — dn4f/dP, are of the same order that those re- 
sulted from the study of magnetovolume effect in SmBg 
(2 Mbar~^ [13]) and from the measurements of resonant 
inelastic X-ray emission in YbAl2 under pressure (~ 5 
Mbar-i dl). 



B. Temperature dependence 

In a simple empirical model which includes interconfig- 
uration fluctuations between /""""^ and /" levels [l^, the 
contribution of the 4/° (J=0) and 4/^ (J=5/2) states of 
Ce to magnetic susceptibility is given by 



X4f(T) = 7VA/i'"4f(T)/3fc(r + Tf). 



(2) 
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Here N is the Avogadro number, fj, effective magnetic mo- 
ment of the 4/^ state, Tf the characteristic temperature 
(valence fluctuation temperature, or Kondo temperature, 
or heavy- fermion bandwidth). It should be mentioned 
that a quantitative analysis of the xaC^) dependence 
using Eq. ^ requires the complete data on n4{{T) (and 
probably on T{{T) as well) which are actually unavail- 
able. Furthermore, to separate the X4f (^) term from the 
experimental data on x{T) one needs to know a back- 
ground contribution xo^ which generally can not be ne- 
glected. A simplified analysis of the experimental data 
can be performed assuming nn, T[ and xo to be tem- 
perature independent. Then the magnetic susceptibility 
obeys a modified Curie- Weiss law. 



X{T) = Xo + X4f{T) = Xo + 



C 



(T-e) 



(3) 



with C = iV/i^n4f/3fc and 8 = — Tf. For the representa- 
tive Ce(Nio.5Cuo.5)5 alloy, the best fit of Eq. ([3]) to the 
experimental data at T > 50 K (Fig. 4(a)) [4I is obtained 
with Xo = 0.6 X 10^"^ emu/mole, C — 0.48 K-emu/mole 
and Q = —79 K. It should be pointed out that the es- 
timate Uif = 0.6, resulted from C, is in a reasonable 
agreement with the value of 0.8 that follows from the 
data in Fig. 2(b) for x ~ 0.5. 

As is evident from Eqs. ([2]) and ([3]), the pressure effect 
on the 4/ susceptibility is governed by changes of n^t and 
Tf with pressure, as 



dlnx4f(T) dlnC 



1 



dP 



X 



dTf 



dP (T + Tf ) dP 
d lnn4f 



dP 



X4f(T) dTf 
C ^ dP' 



(4) 



being a linear function of 1/(T + Tf) or X4f(T). The 
data on dlnx4f/dP for Ce(Nio.5Cuo.5)5 alloy were de- 
rived from the measured effect, dlnx/dP, in the frame- 
work of Eq. ^ by using a value B ^ 1 Mbar"-'^ as a 
rough estimate for the pressure dependence of the back- 
ground [13], which is assumed to originate from 3d{4:d) 
itinerant electrons. The obtained values are plotted in 
Fig. 4(b) as a function of X4f(T). Its linear approxima- 
tion in accordance with Eq. (|4]) gives 



dlnC dlnn4f „„ , _i 

— = , = -3.2 ± 0.7 Mbar \ 

dP dP 

dTf 



dP 



= 1650 ± 250 K • Mbar" 



(5) 



The resultant value dri4f/dP ~ —2.5 ± 0.5 Mbar~ is in 
line with the value dn4f — —2.0 ± 0.3 Mbar"'^ obtained 
previously for x = 0.5 from analysis of the concentra- 
tion dependence of the pressure effect within Eq. ([1]). 
From the pressure dependence of Tf the corresponding 
Griineisen parameter, fi, is estimated as 



dlnT 



dlnV 



B 



dlnTf 
dP 



31 ±5 



(6) 
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FIG. 4: Temperature dependence of the magnetic susceptibil- 
ity X (a) and pressure derivative dlnx4f/dP plotted against 
X4f (b) for Ce(Nio.5Cuo.5)5 alloy. 



using the experimental bulk modulus P = 1.5 Mbar [2l| . 
The Anderson impurity model provides the Kondo tem- 
perature and its pressure derivative to be described in 
terms of n4f f23 |: 



Tk oc 



1 - ?^4f ) d IuTk 



1 



n4f 



dP 



1 - n4f 



d lnn4f 
dP ■ 



(7) 



Then, assuming Tf oc Tk and using in Eq. ([7]) the values 
dlnn4f/dP = -3.2±0.7 Mbar^^ and n4f = 0.8 evaluated 
above for the alloy with x ~ 0.5, one obtains 



^1 f — rij. 



dlnTK 
' d\nV 



24 ±5 



(8) 



in reasonable agreement with the direct estimate 

For Ce(Nio.4Cuo.6)5 alloy, the analogous analysis in the 
framework of Eq. Q and Eq. Q yields the following 
Curie- Weiss parameters: C ~ 0.806 K-emu/mole, xo ^ 
0, Tf = — O = 26 K, and their pressure derivatives: 

dlnC ^ dlnn4f ^ f 

dP dP 
dTf 

— ^ = 620 ± 100 K • Mbar"^ . 
dP 

The latter results in the Griineisen parameter Qf = 35±6, 
assuming the bulk modulus value B = 1.5 Mbar, as in the 
Ce(Nio.5Cuo.5)5 alloy. The similar estimate follows from 
Eqs. ([7]) and © with n4f = 0.93 derived from the data in 
Fig. 2(b). The reasonable description of the Griineisen 
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parameter for alloys with x = 0.5 and 0.6 with Anderson 
model nil allows to consider the Cu-rich alloys studied 
in the present work as the nonmagnetic Kondo lattices. 

V. CONCLUSIONS 

The pressure effect on magnetic susceptibility of 
Ce(Nii_2;Cu2;)5 alloys has been observed for the first 
time. This effect is negative in sign, and also strongly and 
non-monotonously dependent on the Cu content. For the 
reference CeNis compound, the pressure effect value is 
successfully described within LSDA approximation, us- 
ing the modified full potential relativistic FP- LMTO 
method. For Ce(Nii_a;Cuj:)5 alloys the effects of pres- 
sure and alloying on the valence state of Ce ion are the 
most pronounced around x ~ 0.4, which corresponds to 
the half-integer valence v ~ 3.5. In other words, the frac- 
tional occupation nn ^ 0.5 with the nearly degenerate 
/° and configurations of electronic states is favorable 



for the valence instability. It is also found that the main 
contributions to the pressure effect on magnetic suscep- 
tibility for the Cu-rich alloys are i) the decrease of the 
effective Curie constant and ii) the increase of the char- 
acteristic temperature Tf. The latter exhibits a large 
and positive value of the Griineisen parameter, which 
can be apparently described within the Anderson impu- 
rity model. Both of these contributions have their origin 
in the change of the Ce valence state caused by depopula- 
tion of the / state under pressure due to its shift relative 
to the Fermi energy. 
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